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Abstract. We study the possibility of the helium decoupling in the stellar wind of a Ori E. To obtain reliable wind parameters 
for this star we first calculate an NLTE wind model and derive wind mass-loss rate and terminal velocity. Using corresponding 
force multipliers we study the possibility of helium decoupling. We find that helium decoupling is not possible for realistic 
values of helium charge (calculated from NLTE wind models). Helium decoupling seems only possible for a very low helium 
charge. The reason for this behaviour is the strong coupling between helium and hydrogen. We also find that frictional heating 
becomes important in the outer parts of the wind of a Ori E due to the collisions between some heavier elements and the passive 
components - hydrogen and helium. For a metallicity ten times lower than the solar one both hydrogen and helium decouple 
from the metals and may fall back onto the stellar surface. However, this does not explain the observed chemical peculiarity 
since both these components decouple together from the absorbing ions. Although we do not include the effects of the magnetic 
field into our models, we argue that the presence of a magnetic field will likely not significantly modify the derived results 
because in such case model equations describe the motion parallel to the magnetic field. 
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1. Introduction 

The radiative force plays an important role in the atmo- 
spheres of early-type stars. In luminous early-type stars it may 
be greater than the gravitational force. In such a case, the 
outer stellar layers can hardly be in hydrostatic equilibrium. 
Consequently, the radiatively driven wind is blowing from 
these stars (see, e.g., Kudritzki & Puis 2000] Owocki l2"0Tj4l or 
Krticka & Kubat 2006 for a review). For many cooler early- 
type stars the radiative force is not able to expel the mate- 
rial into the circumstellar environment, however it is able to 
cause elemental drift and, consequently, abundance stratifica- 
tion (e.g. Vauclair 2003, Michaud 2005) and becomes mani- 
fested as chemical peculiarity. 

Hot stars with radiatively driven stellar winds and hot 
chemically peculiar stars seem to form two diverse groups of 
stars. While the flow in the atmospheres of stars with strong 
winds would eliminate any possible abundance stratification 
and chemical peculiarity, the atmospheres of chemically pecu- 
liar stars are supposed to be very quiet (likely due to magnetic 



fields) to enable elemental diffusion (e.g. Vauclair fl975 



However, there exists a small group of hot stars for which 
both effects, stellar wind and chemical peculiarity, seem to be 
important. The sparsely populated group of helium rich stars 
probably exhibits both effects. In some cases, chemical pecu- 
liarity can be even used as a test of wind existence (Landstreet, 
Dolez, & Vauclair fTWSI Dworetsky & Budaj|3000). Some of 
the stars with parameters corresponding to hot chemically pe- 
culiar stars may have also purely metallic winds, as proposed 
bv Babel (fT996b . 

One of the most enigmatic stars of the group of helium rich 
stars is the star denoted as <x Ori E or HD 37479. To our knowl- 
edge, this star is the only known Bp star with hydrogen lines in 
emission. Its photometrical and emission line profile variabil- 
ity can be explained assuming the circumstellar cloud model 
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(Smith & Groote l20Qll Townsend et al. 120051 . Similarly to the 
star HD 37776 (Mikulasek et al. I2006l > this star may exhibit 
rotational braking probably due to the magnetically controlled 
stellar wind (Oksala & Townsend 2006 1. 

Stellar winds of hot stars are accelerated mainly due to the 
light scattering in the lines of elements like carbon, nitrogen, 
oxygen or iron. Because these heavier elements have much 
lower density than the bulk of the flow (composed of hydro- 
gen and helium), hot star winds have a multicomponent nature 
(Castor, Abbott, & Klein [l976l Springmann & Pauldrach[T332| 
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Cure fTff^l Babel [T995l Krticka & Kubat l2UU0ll2Wn hereafter 
KKII). Momentum is transferred from low-density heavier el- 
ements to the high-density component composed of hydrogen 
and helium due to Coulomb collisions of charged particles. The 
corresponding frictional force depends on the velocity differ- 
ence between wind components. For very low velocity differ- 
ences, lower than the mean thermal velocity (basically for stars 
with dense winds), the flow is well coupled. For higher veloc- 
ity differences (basically for stars with low density winds) the 
Coulomb collisions are not able to efficiently transfer momen- 
tum from heavier ions to hydrogen and helium, and frictional 
heating becomes important (KKII) or wind components may 
even decouple (Owocki & Puls l20U2l Krticka & Kubat 12002b . 

The chemical peculiarity of helium strong stars is in the 
literature explained either by radiative diffusion moderated by 
the stellar wind (Vauclair [T975l Michaud et al. 1987 1 or by he- 
lium decoupling in the stellar wind and its consecutive fall back 
(Hunger & Groote 1999 ). To test whether the wind decoupling 
model is able to provide a reliable explanation of helium chem- 
ical peculiarity we decided to calculate multicomponent wind 
models suitable for the Bp star cr Ori E. 

2. NLTE wind model 

2. 1 . Basic assumptions 

Since there is no accurate determination of cr Ori E wind pa- 
rameters from observation, we decided first to calculate wind 
mass-loss rate and terminal velocity using our NLTE wind 
models (Krticka & Kubat l2004l 

These models make it possible to solve (generally multi- 
component) hydrodynamical equations for radiatively driven 
winds and enable prediction of mass-loss rates and terminal 
velocities. The radiative force and the radiative heating/cooling 
term are calculated with NLTE level populations obtained us- 
ing atomic data based on an extended set of the TLUSTY 
input files (Hubeny IT9881 Hubeny & Lanz fT992l Hubeny & 



Lanz 1995 Lanz & Hubeny 2003). The data are based on the 
Opacity Project (Seaton Luo & Pradhan [T989I Sawey 

& Berrington, Seaton et al. [T9"9"2l Butler et al. [T9"9"3l Nahar & 
Pradhan[1993 1, the Iron Project (Hummer et al. 119931 Bautista 
[T996I Nahar & Pradhan fT9%l Zhang fT99TI Bautista & Pradhan 
[T9971 Zhang & Pradhan [T997I Chen & Pradhan fT999l and the 
VALD database (Piskunov et al. fl995l Kupka et al. fT999l . The 
most important simplification of our code compared to other 
models available in the literature ( Vink et al. 119991 Pauldrach et 
al. 2001 Grafener & Hamann 2005 ) is the simplified treatment 
of the radiative transfer equation. We use Sobolev approxima- 
tion for the solution of the radiative transfer equation in lines 
(Sobolev fT947l Rybicki & Hummer fT978l and we neglect UV 
line blocking as well as line overlaps. 

For the calculation of wind models we assume spherical 
symmetry. Since it is well-known that cr Ori E has a magnetic 
field (Landstreet & Borra 1978 1 that generally breaks down the 
assumption of spherical symmetry, our models, at least approx- 
imately, describe the flow along individual field lines. Note that 
in order to obtain a more realistic description of the multicom- 
ponent wind in presence of a stellar magnetic field, it would 




Fig. 1. The radial variation of helium ionization fractions in 
the wind calculated using NLTE model. He n is the dominant 
ionization stage throughout the wind. Close to the star (above 
the photosphere) the wind density decreases, consequently the 
wind becomes more ionized and the ionization fraction of He i 
decreases with the increasing radius, while the fraction of He in 
increases with the radius. In the outermost regions the wind 
becomes slightly less ionized again due to the decrease of the 
mean intensity. 



be necessary to solve the hydrodynamical equations along field 
lines with inclusion of proper projections of individual forces 
acting on the wind or to use the magnetohydrodynamical simu- 
lations. In the case of a one-component wind model the former 
was done in a simplified form by Friend & MacGregor ( 1984 1, 
MacGregor & Friend ( 1987 1, and more recently the latter using 
magnetohydrodynamical simulations by ud-Doula & Owocki 
<2002|). 

Stellar parameters appropriate for cr Ori E (Hunger et 
al. 0^89) are given in Table Stellar radiative flux (lower 
boundary of the wind solution) of cr Ori E is calculated us- 
ing the spherically symmetric NLTE model atmosphere code 
described in Kubat (2003). 

Since we intend to study self-initiation of the chemical pe- 
culiarity of cr Ori E, we assume solar helium abundance. If not 
explicitly stated, we assume solar chemical composition also 
for heavier elements. 



Table 1. Adopted stellar (Hunger et al. [1989) and calculated 
wind parameters of cr Ori E. 



effective temperature r e fr 


22 500 K 


stellar mass M 


8.9 M G 


stellar radius R, 


5.3 R Q 


logg[CGS] 


3.94 


mass-loss rate M 


2.4 x lO^Moyr 1 


terminal velocity 


1460 km s" 1 
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2.2. Calculated NLTE model 

Using our NLTE wind model, in which we neglected the multi- 
component effects, we calculated wind parameters for cr Ori E. 
The predicted wind mass-loss rate M = 2.4 x 1CT 9 M Q yr _1 is in 
agreement with order of magnitude estimate inferred from ob- 
servations (see the discussion in Groote & Hunger 1997 1. The 
derived terminal velocity Voc = 1460 km s _I implies the ratio of 
terminal velocity to escape velocity Voo/v esc = 1 .8 (note that the 
radiative force due to light scattering on free electrons has only 
a small influence here since the Eddington parameter F « 0.02 
is low). The derived ratio of terminal velocity to escape veloc- 
ity is consistent with that for stars close to the bistability jump 
(see Pauldrach & Puls[l990 Lamers et al. 119951 Crowther et 
al. l2006l . 

Since we want to study helium decoupling that is charge- 
dependent (as follows from the frictional term due to the 
Coulomb collisions - see Eq. (|6j or Eq. 40 in KKII) we have 
to study the calculated variation of He ionization in the wind. 
From Fig.^we conclude that helium is mostly singly ionized 
in the stellar wind with only a small contribution of doubly ion- 
ized helium. The fraction of neutral helium in the wind is very 
small, less than 10~ 4 (for r > 1.01/?.). 



0.08 




1 1.5 2 2.5 3 3.5 4 4.5 



Fig. 2. Calculated relative velocity differences (see Eq. flj) be- 
tween individual heavier elements and the passive component. 
Velocity differences for elements with low x/, p values are not 
plotted in this graph. Note that here the effects of frictional 
heating are neglected. 



2.3. Multicomponent effects 



Before studying helium decoupling we shall discuss another 
aspect of a wind model for cr Ori E, namely, multicomponent 
effects due to the friction of individual heavier elements with 
hydrogen and helium. 

Since the radiative acceleration acting on individual heav- 
ier elements is different, each element moves with slightly dif- 
ferent velocity. The dimensionless velocity difference between 
given individual heavier element h and hydrogen and helium 



Q7,p 



is roughly equal to (Krticka 2006 1 



Xkp ~ 8h 



3kT 



«p 8V^^plnA 



(1) 



(2) 



We use the subscript h to denote individual heavier elements 
- i.e. h stands for N, Si, S, etc. For this purpose we do not 
use the subscript i that denotes all heavier elements described 
together. The subscript p denotes the passive component (hy- 
drogen and helium). Here v p and v/, are the radial velocities of 
the components, m p and to/, are their atomic masses, q/, and q p 
their charges, g™ d is the radiative acceleration acting on given 
element h, n v is the number density of hydrogen and helium, T 
is the temperature (assuming that the temperatures of all wind 
components are the same), In A is the Coulomb logarithm, and 
finally, 

2kT [mp + to/,) 



nipirih 



(3) 



Using Eq. (|2ji we calculated the approximate velocity dif- 
ferences between hydrogen and helium and individual heavier 
ions (see Fig. [3}. Apparently, the relative velocity differences 



between silicon or nitrogen and the hydrogen-helium compo- 
nent are so high that frictional heating may be enabled. To test 
this possibility we calculated four-component wind models. 
Despite the fact that the relative velocity difference between 
silicon and the passive component is slightly higher than that 
of nitrogen, we selected the nitrogen component for this test be- 
cause nitrogen has higher abundance and, consequently, heat- 
ing effects caused by friction are expected to be higher for this 
element. Resulting four-component wind model (with wind 
components hydrogen plus helium, free electrons, nitrogen and 
the remaining heavier ions) is given in Fig. [5] Apparently, the 
stellar wind temperature significantly increases in the outer 
parts due to frictional heating. However, no decoupling occurs 
in this case (since xn p is too low to enable the decoupling). The 
wind terminal velocity significantly increases compared to the 
case when the frictional heating is neglected. This is caused by 
the change in the ionization and excitation equilibrium due to 
a higher wind temperature. Since this effect occurs for veloc- 
ities higher than that at the critical point, the mass-loss rate is 
not essentially modified in four-component models. A similar 
effect, but to a lesser extent, was already reported by Krticka 
d2006b . 

Because the effects of frictional heating between heavier 
elements (in our case nitrogen) and the passive component are 
important basically only for velocities higher than the escape 
velocity (which is v esc = 790 kms -1 at the stellar surface), they 
do not significantly modify the flow properties close to the stel- 
lar surface. Since this region is important for the study of the 
helium decoupling and its reaccretion, we neglect the effect of 
friction caused by the distinct properties of all individual heav- 
ier elements. However, heating due to friction between hydro- 
gen and helium ions is properly taken into account. 
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Fig. 3. Four component NLTE wind model with components 
hydrogen and helium (p), electrons, nitrogen (N) and remain- 
ing heavier ions. Upper panel: Radial velocity of the passive 
component p and nitrogen. The velocities of wind components 
are nearly indistinguishable in the plot. Middle panel: Plot of 
dimensionless velocity difference between the passive compo- 
nent and nitrogen x Np (Eq.[T). The velocity differences between 
wind components are relatively high, enabling frictional heat- 
ing (as shown in the plot of the temperatures of individual wind 
components, lower panel). 



3. The possibility of helium decoupling 

3.1. Model assumptions 

To test the possibility of helium decoupling we calculated four- 
component wind models with the following four explicitly in- 
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Fig. 4. Calculated velocity structure of a four-component (H, 
He, heavier ions, electrons) wind model for a realistic value 
of helium charge z a — 1 . Hydrogen and helium velocities are 
nearly identical and helium does not decouple. 



eluded wind components: hydrogen, helium, heavier ions and 
free electrons. For this purpose we did not use our NLTE mod- 
els, but simpler (and faster) models with the radiative force 
in the CAK approximation (Castor, Abbott & Klein 1975J see 
KKII for the description of these models). The parameterized 
description of the radiative force in these models enables us to 
vary both wind parameters and wind ionization state indepen- 
dently to discuss the parameter space for which helium decou- 
pling seems possible. 

To model mass and charge of heavier ions we selected pa- 
rameters corresponding to oxygen (because elements like C, 
N, O, Si, S are important for radiative driving). The selection 
of oxygen does not have a significant impact on derived results. 
The CAK force parameters 



£ = 0.38, 



0.50, 



5 = 0.03 



(4) 



(CAK, Abbott [T982I Puis et al. EOOOl l are chosen to obtain the 
same mass-loss rate and a similar velocity profile as derived 
from our NLTE models. Note that the terminal velocity calcu- 
lated using these parameters is slightly higher than calculated 
from NLTE models. However, this does not significantly influ- 
ence our analysis. 

3.2. Realistic helium charge 

We calculated four-component wind models with the realistic 
value of the helium charge taken from previously discussed 
NLTE model. Since helium is mostly singly ionized in the wind 
region (see Fig.^, we simply assume z a = 1, where z a is the 
relative charge in units of the elementary charge e, q a = ez a . 

The velocity profile of a calculated four-component wind 
model is shown in Fig.|4] Evidently, the radial velocities of all 
wind components are nearly the same and the wind compo- 
nents are well coupled. Consequently, helium decoupling can- 
not be found in this case. 
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Fig. 5. The velocity structure of a calculated four-component 
wind model for a very low probably non-realistic value of he- 
lium charge z a = 0.01. Helium is accelerated at the wind base, 
however due its low charge it subsequently decouples and may 
fall back. 

3.3. Low helium charge 

For very low helium charge, significantly lower than that calcu- 
lated from NLTE equations, helium decoupling is possible. As 
an example we plot in Fig.|5]a comparable wind model but for 
helium charge z a - 0.01. Apparently, the wind density at the 
wind base is high enough to accelerate helium from the stel- 
lar surface to supersonic velocities. However, due to its very 
low charge helium decouples from the mean flow at veloci- 
ties lower than the escape velocity. Consequently, helium re- 
mains bound to the star and may fall back to the stellar surface. 
Hydrodynamical simulations are necessary to study this pro- 
cess in detail (cf. Porter & Skouza lT^I Votruba et al. 2006 1. 



3.4. Helium coupling 

The main result of this paper - that for realistic values of its 
charge helium does not decouple from the wind flow for pa- 
rameters suitable for a star like <x Ori E - seems to be in contra- 
diction to the result of Hunger & Groote ( 1999 ) who argued on 
the basis of the work of Babel ( 1995 1996 1 for the decoupling 
hypothesis that explains observational results. To find the rea- 
son for this difference we study the hydrodynamical equations 
analytically. The equation of motion for helium (neglecting the 
radiative force acting on helium, the gas-pressure term and the 
charge separation electric field) has the form of (Burgers 1969 1 

Va-T- = ~8 + — > . K aa Gch(Xaa)] T, (5) 

where v a is the helium radial velocity, p a is the helium mass 
density, Gch(x) is the Chandrasekhar function, and g is the 
gravitational acceleration, g = GM/r 2 . The summation of the 
frictional force in Eq. Q takes into account friction due to hy- 
drogen "H", free electrons "e" and all heavier ions "i" ( the lat- 
ter are described as one component). The frictional coefficient 
K ab due to the Coulomb collisions between any component a 



W a q 2 b , . 

K ah = n a n b —— — In A, 

kT ah 



(6) 



where n a and tif, are the number densities of individual com- 
ponents and q a , q b are charges of these components. The mean 
temperature is 

m a T b + m b T a 

i a b = • (') 

m a + m b 

In the case of the helium decoupling, helium has to decou- 
ple from both hydrogen and heavier ions. The maximum pos- 
sible value of the Chandrasekhar function is in both cases the 
same, however, the value of the frictional coefficient is differ- 
ent, K a u » K al since «h » n\ (assuming an ionized medium). 
This means that hydrogen-helium collisions are extremely im- 
portant for the dynamics of helium component. Consequently, 
hydrogen may be (depending on its ionization state) much more 
important for the helium acceleration than the heavier ions. The 
neglected strong coupling between hydrogen and helium is the 
main reason why Hunger & Groote ( 1999 1 concluded that the 
helium decoupling from flow of hydrogen and passive plasma 
is possible (even for realistic parameters of wind plasma - like 
helium ionization state). 

However, the situation may differ from the extreme case de- 
scribed above. The contributions of hydrogen and heavier ions 
to helium acceleration may be equally important due to differ- 
ent values of the relative velocity differences x a n and x a \ that 
are arguments of the Chandrasekhar function (see Appendix 
lAl for a more detailed approximate description of the flow). 
However, in the case when the velocity difference between he- 
lium and heavier ions is large, x a \>\, the frictional force due 
to hydrogen dominates. Note that after the decoupling of he- 
lium and hydrogen (in the particular case displayed in Fig. [5] 
for r//?» > 1.06), helium still remains coupled to electrons. 
This is due to the fact that a ae » a a n- The coupling of he- 
lium to electrons, however, modifies the helium velocity only 
slightly. 

The only possibility how to enable helium fall back onto 
the stellar surface is to break the strong coupling between he- 
lium and hydrogen while keeping coupling between hydrogen 
and heavier ions. The only way how to secure this in our self- 
consistent models is to lower the helium charge. This is why 
for very low helium charge helium decoupling is possible. 

We may derive an approximate condition for the helium 
decoupling in the stellar wind. If helium decouples from the 
stellar wind and subsequently falls back, there is a point where 
the helium velocity gradient is zero (see Fig. |5jl- The helium 
momentum equation Eq. has in this point the form of 



KHaGch(XHa) 

Pa 



(8) 



(neglecting friction with heavier ions and electrons). The 
Chandrasekhar function is limited by its maximum value 
Gch(xna) ^ G max . In the case of decoupling the frictional 
acceleration on the right hand side of Eq. (|8} is lower than 
the gravitational acceleration (inequality in © occurs). With 
g = GM/r 2 , the hydrogen mass-fraction X as 0.71, and using 
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Eq. (|6) the condition for helium decoupling for the radial ve- 
locity v is 

kGM . Mz- 



e 4 XlnAG n 



vT 



(9) 



where e is the elementary charge. Note that the wind tempera- 
ture T (or, explicitly, the mean temperature T Ho , see Eq.Q) is 
in our case usually influenced by the frictional heating. Using 
scaled quantities, (|9) can be rewritten as 



M 



[10-™ Mo/year; 



zl (lMo/do^doOkms- 1 )' (10) 



A similar condition was obtained by Owocki & Puis (2002 
Eq. (23) therein), but for coupling between hydrogen and heav- 
ier ions. The inequality (II 01 gives condition for the helium 
charge that enables helium decoupling. For example, for values 
taken from the model displayed in Fig.|5]we obtain the relative 
helium charge necessary for decoupling z a ~ 0.01, consistent 
with our previous detailed calculations. 

3.5. The possibility of the decoupling of neutral helium 

Our conclusion that helium decoupling is unlikely in the case 
of realistic helium charge is based on the assumption that all 
helium ionization stages can be treated as a single wind com- 
ponent. This physically means that the processes of ionization 
and recombination are so rapid that the typical time of re- 
combination or ionization is much shorter than the character- 
istic time of change in flow properties (i.e. velocity and den- 
sity). While the violation of this condition would not affect 
the results in the case of doubly ionized helium (because it 
has higher charge than the singly ionized helium and, conse- 
quently, is more tightly coupled to hydrogen, see Eq. [6ji, the 
violation of this condition in the case of neutral helium may 
offer an interesting possibility of a helium decoupling mecha- 
nism. The reason for this is that the cross-section of collisions 
of protons with neutral helium is much lower than the effective 
cross-section for collisions between charged particles (Krstic 

& shuitz figggi 

To test this, we compare the characteristic time of change 
in flow properties (roughly corresponding to the time during 
which the velocity increases to its actual value) 



T flow - T" 



with the typical time of ionization of He i 

1 



(ID 



(12) 



"He HHe II 

where PnenHeii is the ionization rate of Hei (Mihalas 1978 1. 
This value was taken from our NLTE models. Since excited 
levels of neutral helium have very low occupation numbers, we 
neglected their contribution to ^Hei^Hen- 

Comparison of characteristic times rg ow and Ti on i z is given 
in Fig. [6] The characteristic time ra ow first slightly decreases 
with radius and then increases again. This marks the region 
of fast wind acceleration. Also the Hei ionization time in- 
creases with radius. However, the characteristic time Tfl OW is 
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Fig. 6. Comparison of the characteristic time of change in flow 
properties Tfl OW and the characteristic time of Hei ionization 

Tioniz 

much larger than the Hei ionization time Ti on i z . Consequently, 
helium may be treated as one component and no decoupling 
between different helium ionization stages occurs. 

4. Models with lower metallicity 

Since the metallicity of individual surface elements of cr Ori E 
varies with location on the stellar surface, we studied also 
the influence of different metallicity on the possibility of he- 
lium decoupling. In particular, we selected the metallicity to be 
Z/Z Q = 0.1, as derived by Smith & Groote (USD for cr Ori E. 

4.1. One-component NLTE wind model 

Our NLTE wind model for metallicity Z/Z = 0.1 (not shown 
here) predicts a very low mass-loss rate of M = 7.4 x 
10~ 11 M yr _1 . This is much lower than that can be inferred 
from simple relation M ~ Z 67 derived by Krticka (2006 1 
for SMC O stars. This indicates that for lower wind densities 
the dependence of mass-loss rate on metallicity is more pro- 
nounced. 

Steeper metallicity dependence of mass-loss rates for small 
metallicities may seem unexpected. To explain this behaviour, 
let us first discuss the contribution of lines of different ele- 
ments to the radiative force. This is connected with the value 
of the line optical depth t which is in the Sobolev approxima- 



tion given by (e.g. Abbott 1982 , Eq. 1, neglecting the finite-disk 
correction) 



m e Vij \g 



si 



dv 



dr 



(13) 



where n,, ttj, gt, gj are the number densities and statistical 
weights of individual states i, j that give rise to the studied line, 
gi fij is the line oscillator strength and, vy is the line frequency. 
The optically thick lines (t > 1) are important for the radiative 
acceleration. Note that the radiative force due to the optically 
thick lines does not depend on number densities of involved 
atomic states. 
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Very close to the star, the Sobolev optical depth Jl 3i is 
large for many lines both of iron group elements and of lighter 
elements (like carbon, nitrogen, etc.) mainly due to the very 
large wind density. Since iron group elements have effectively 
larger number of lines than lighter elements, the former dom- 
inate the radiative force close to the star. Farther out in the 
wind the line optical depth Jl 3i is lower (mainly due to de- 
creasing wind density). Since iron group elements have lower 
abundances than carbon and nitrogen (i.e. effectivelly lower n,), 
many e.g. iron lines become optically thin (t < 1) and these 
lines do not significantly contribute to the radiative force for 
higher wind velocities. Consequently, the relative contribution 
of iron group elements to the radiative force decreases with in- 
creasing wind velocity and lighter elements become more im- 
portant for the radiative acceleration (see also Vink et al. 2001 1. 
This behaviour can be also interpreted as the result of the fact 
that lighter elements and iron group elements dominate dif- 
ferent parts of the line-strength distribution function (Puis et 
al. 2000). In our case (even for the solar metallicity model) the 
iron lines become inefficient already close to the critical point 
due to a very low wind density. 

However, for low metallicities and densities also many lines 
of lighter elements that are for higher metallicities optically 
thick, become optically thin. This leads to a significant de- 
crease of the mass-loss rate. In our specific case there is only 
one optically thick line ( the C m resonance line at A = 977 A) 
that significantly contributes to the radiative force at the criti- 
cal point. From the contribution of this optically thick line we 
can calculate the value of the a parameter (cf. Puis et al. 2000 
Eq. 26) as a = 0.46. Using Eq. (87) of Puis et al. (HHP, we 
can roughly estimate the dependence of the mass-loss rate on 
the metallicity as M ~ Z 1 1 . This is more consistent with the 
decrease of the mass-loss rate with metallicity obtained in our 
case than the dependence M ~ Z 67 derived by Krticka (2006) 
for the denser wind of SMC stars. Note however that the value 
of a may itself depend on the metallicity. In the terms of the 
line-strength distribution functions (Puis et al. 2000) the steep- 
ening of the dependence of mass-loss rate on metallicity may 
be explained as a consequence of the steepening of the line- 
strength distribution function for high line-strengths. Our re- 
sult supports the findings of Kudritzki ( 2002 see Fig. 9 therein) 
who also found more pronounced dependence of mass-loss rate 
on metallicity for a very low-metallicity environment. 

The calculated wind model is indeed very close to the wind 
existence limit. As discussed above, iron lines become ineffi- 
cient for wind driving and the stellar wind itself is accelerated 
mainly due to carbon and nitrogen lines. CAK critical point is 
located closer to the stellar surface as a result of the decrease of 
a number of optically thick lines (cf. Babel l 19951 Vink l2006> . 

For the calculation of NLTE models we neglected the multi- 
component effects, however from the values of relative velocity 
differences it is clear, that the stellar wind would decouple 
relatively close to the stellar surface. We have tested this pos- 
sibility in the three-component model where the metallic ions 
are described as one component and concluded that already in 
this model the decoupling occurs close to the stellar surface. 

The wind terminal velocity v M = 1 120kms _1 is lower than 
in the case of a solar-metallicity model (cf. Fig.©. The lower- 
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Fig. 7. The velocity structure of calculated four-component 
wind model with low metallicity Z/Z Q = 0.1 and a realistic 
value of helium charge z a — 1 ■ Helium velocity profile is nearly 
indistinguishable from the hydrogen one. Hydrogen and helium 
are accelerated at the wind base, however due to low wind den- 
sity they decouple from heavier ions and may fall back. 



ing of the terminal velocity with respect to the solar metallicity 
model is an example for the dependence of the terminal veloc- 
ity on metallicity. 

The most important parameter that influences the possibil- 
ity of helium decoupling is the helium charge. However, due to 
the lower wind density, helium is even more ionized than in the 
case of the solar metallicity model. 

4.2. Four-component CAK model 

We have calculated low-metallicity four-component wind mod- 
els with helium as the fourth component again with the sim- 
pler models using the CAK radiative force. Although the he- 
lium ionization is slightly higher than that in the solar metal- 
licity model, we assumed z a = 1. To obtain the lower mass- 
loss rate inferred from NLTE calculations we used k = 0.23. 
The other force multipliers are the same as in Eq. ©. It may 
seem that the adopted value for k was not adequatelly low- 
ered compared to the value adopted for solar metallicity model 
(where k = 0.38 was used), because k scales with metallicity 
as k ~ Z l - a (e.g. Puis et al. l2000l . However, the value k = 0.23 
is appropriate for a four-component wind model of a star with 
low metallicity (i.e. Z/Z = 0.1). Because the radiative force in 
the four-component models is calculated using metal density, 
the decrease of wind metallicity causes lower radiative force 
and lower mass-loss rate. Consequently, with decreasing metal- 
licity the value of k in a four-component models varies only 
slightly. 

The results from the calculated four-component wind 
model can be seen in Fig.Q Apparently, the friction between 
metallic ions and the passive component is low and is not able 
to accelerate the passive component to velocities higher than 
the escape velocity. Consequently, the wind components de- 
couple and the passive component may start to fall back to the 
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stellar surface. However, hydrogen and helium are still coupled 
together. 

From the first point of view it could be expected that 
there would be a very high velocity gradient of heavier ele- 
ments downstream from the point where the decoupling occurs. 
However, this is not the case, the velocity increase of heavier 
elements is only moderate after the decoupling. The reason for 
this behaviour is the same as that which causes non-decoupling 
described for the two-component flow by Krticka & Kubat 
( 2000 see the arguments in the Chapter 3.2 therein). According 
to the momentum equation of heavier ions, the radiative force 
is approximately equal to the frictional force. Since there is 
a deficit of the frictional force in the region where decou- 
pling starts to be effective, the radiative force is lower there. 
Consequently, due to the dependence of the radiative force on 
the velocity gradient in the Sobolev approximation, the veloc- 
ity gradient of heavier ions decreases. Because the velocity of 
the passive component (hydrogen and helium) is lower than the 
escape velocity, this component decelerates and may fall back 
onto the stellar surface (Porter & Skouza 1999 1. The wind in 
this region is likely unstable for the perturbations of a multi- 
component flow (Owocki & Puls l2002l Krticka & Kubat l2002> . 

The derived terminal velocity for the absorbing ions is very 
low in this case, roughly 130kms _1 . Consequently, due to low 
wind velocities there would be basically no X-rays emitted in 
this case by the possible wind shocks (provided that the ab- 
sorbing ion component does not accelerate as a result of other 
instabilities). This model corresponds to models discussed by 
Groote & Schmitt ( 2004 1, who argued that the possible expla- 
nation of missing periodic X-rays comming from cr Ori E may 
be connected with too low wind velocities. 

5. Discussion 

Although we have used the best models available to us, there 
are many simplifications that may slightly modify the results 
obtained (see also KKII, Krticka 2006 1. 

Mass-loss rate: The mass-loss rate of cr Ori E may be actually 
lower than that obtained from NLTE wind models, similarly 
as deduced from observations of O stars with low luminosities 
(e.g. Bouret et al. 120031 Martins et al. I2004t . To obtain he- 
lium decoupling for cr Ori E due to a low mass-loss rate, mass- 
loss rates of order 10~ 13 M year -1 are necessary. However, in 
such a case also hydrogen decouples from the flow (e.g. KKII). 
Therefore, we conclude that the change of the wind mass-loss 
rate likely does not influence our basic results. 

Helium-rich stars could also eventually have a purely 
metallic wind, as proposed by Babel ( 1996). However, it is not 
clear how this could lead to helium chemical peculiarity. 

Clumping: Winds of OB stars may be clumped (see Martins 
et al. (20051 or Fullerton et al. (2006 1 for a recent observa- 
tional support of the existence of wind clumping in O stars). 
It is not clear how this effect influences the multicomponent 
flow. Whereas the physical state of the interclump medium 
may significantly change due to the processes connected with 



the decoupling, it is possible that this effect is not very sig- 
nificant inside the clumps. This is caused by the fact that the 
frictional force depends on the wind density as p 2 , but on the 
other hand the mass-loss rate derived from observations may 
be overestimated due to the clumping. Possibly, both effects 
may mutually cancel and the net effect of wind clumping may 
be marginal. This picture however significantly depends on the 
physical properties of clumps and of surrounding media. 

Magnetic field: We have also neglected the influence of the 
magnetic field although it is known that cr Ori E has a strong 
one. The self-consistent treatment of this problem is very com- 
plicated and involves the solution of magnetohydrodynamic 
equations (e.g., ud-Doula & Owocki 2002 1. It should be noted 
that for such strong field as that of cr Ori E this has not yet been 
done, even in the one-component approximation, due to ex- 
tremely low Courant time steps. An approximate treatment can 
be done by solving the hydrodynamical equations along cor- 
responding field lines (e.g. Mestel 1968). Since this approach 
applies merely a geometric projection of corresponding terms 
in hydrodynamical equations, it is not likely that this would 
change the possibility of helium decoupling significantly. 

Although MacGregor 1 1988) showed that strong magnetic 
field may significantly influence the wind terminal velocity, 
Owocki & ud-Doula ( 2004 1 concluded that the change of the 
terminal velocity is only moderate. The latter argument sup- 
ports our basic results to be valid also in the case of strong 
magnetic field. 

The magnetic field, however, may also change the frictional 
coefficient due to the particle gyration induced by the magnetic 
field or may also alter the velocity distribution function of in- 
dividual species. The modification of frictional coefficient due 
to the magnetic field could explain the observational finding of 
Groote & Hunger ( 1997 ) that the helium-enrichment in cr Ori E 
is found where the windflow in the non-magnetic case will be 
mostly parallel to the magnetic field lines of this star. Similar 
findings for other helium strong stars were reported by Groote 
(120031 

The importance of the magnetic field for the collision fre- 
quency of helium atoms is given by the value of the product of 
the mean collision time between helium ions and protons r a n 
and gyration frequency of helium ions u>„ (Mestel 2003 page 
18 therein; Braginskii 119631 Imazu[1986l. The mean collision 
time r ff H is given by 

T -i - -1 ^liL d4) 

t?H "3V^W,' ] 

where the frictional coefficient K a u is given by Eq. (|6), and a a u 
is given by Eq. Q. The gyration frequency is 

GJ a = , (15) 

m a c 

where B is the magnetic field intensity. For collision-dominated 
plasmas T a nu> a <K 1 the influence of the magnetic field on 
the momentum transport between individual plasma compo- 
nents can be neglected. In the opposite case T a nu> a » 1 a 



J. Krticka, J. Kubat, and D. Groote: Multicomponent radiatively driven stellar winds IV. 



9 



strong anisotropy of microscopic plasma properties may oc- 
cur (in the directions parallel and perpendicular to the mag- 
netic field). However, the frictional coefficients parallel to the 
magnetic field remain the same as in the case without magnetic 
field. Consequently, the change of frictional coefficients due to 
the magnetic field is not important in the case when magnetic 
energy density dominates over the gas kinetic energy density, 
i.e. in the case when plasma moves along magnetic field lines. 

Application of Eqs.^JEJto our derived wind model has 
shown that for typical magnetic field with intensities of the or- 
der of 10 4 G, which corresponds to the surface magnetic field 
of cr Ori E derived by Landstreet & Borra ( 1978 1, in the outer 
parts of cr Ori E wind the inequality t„hw« s> 1 holds. This 
inequality holds even for lower magnetic field intensities of 
the order of 10 2 G (the lower value of the magnetic field in- 
tensity roughly takes into the account its decrease with radius 
and its latitudinal variations). Since wind flows along magnetic 
field lines due to the strong wind confinement by the magnetic 
field, the above mentioned variation of frictional coefficient 
with magnetic field is likely not very important in our case. 

Ionization: Our results are sensitive to the ionization state of 
helium atoms. However, because derived helium charge is by 
two orders of magnitude higher than the value necessary for 
helium decoupling, it seems that our results are also robust 
with respect of errors in the calculation of statistical equilib- 
rium equations (note that helium is dominated by its first ion- 
ization stage already in the stellar atmosphere). 

Even in the case when doubly ionized helium comprises 
a significant fraction of all helium atoms, some of them re- 
main neutral (Fig. [Q. This means that each helium atom for 
some (however very small) fraction of time is neutral. Recall 
that in the case of cr Ori E some wind material may remain 
trapped in the minima of gravitational and centrifugal poten- 
tial along magnetic field lines (Townsend & Owocki 2005 1. 
The non-zero fraction of time spend by each helium atom in 
the neutral state may then induce drift of helium atom across 
magnetic field lines. Let us at least approximatively derive the 
magnitude of such an effect. Recall that the typical time of 
neutral helium ionization is Tioniz (see Eq. H2\ ) and the typ- 
ical time of neutral helium recombination is introduced simi- 
larly as T re c- During the time Ti on ; z the helium atom is not sup- 
ported by the magnetic field (since it is neutral) and its abso- 
lute value of the radial velocity increases by (neglecting the 
centrifugal force) gT lon i z where g = GM/r 2 . During the time 
T rec the helium atom is supported by the magnetic field and has 
zero radial velocity. Assuming that helium is mostly ionized, 
Tree s> Ti on i z , this means that the average velocity is roughly 
£ T ioni z /( T rec + Tioniz) ~ 8 T ^ on J T Kc- The typical time necessary 
to travel the distance R* is Tag - RIt kc / |GMr? oniz j. For typical 
values for the outer wind parts r rec m 10 5 s and Ti on i z ~ 1 s and 
the typical diffusion time t&b ~ 10 5 years is too large compared 
to the typical time for emptying of the magnetosphere, which 
is of the order of a hundred years (Townsend & Owocki 2005 1. 
Moreover, the drift may occur in outward direction, since the 
total force acting on neutral helium atoms (i.e. the centrifugal 
force plus the gravitational force) may be directed outwards 



and in such a case no fall back will occur. Note that the rota- 
tion of the magnetosphere containing charged particles may be 
roughly approximated as a rigid body rotation, at least near to 
the stellar surface. 

6. Conclusions 

We have studied whether the helium decoupling can occur 
in the stellar wind of cr Ori E and, consequently, we tested 
whether this process may consistently explain enhanced he- 
lium abundance in the atmospheres of Bp stars (as proposed 
by Hunger & Groote IT9"9"9"l . 

For this purpose we first calculated an NLTE wind model 
for cr Ori E and derived its mass-loss rate, terminal velocity 
and helium ionization. These parameters were consequently 
used for the calculation of four-component wind models to test 
whether helium decoupling may occur. As a by-product of our 
study we have found that the frictional heating due to colli- 
sions between individual heavier elements (for example nitro- 
gen) and the passive component (hydrogen and helium) is im- 
portant in the outer wind regions of cr Ori E. 

Our four-component wind models showed that for realis- 
tic values of helium charge the decoupling of helium from the 
stellar wind of cr Ori E is unlikely. Helium decoupling is possi- 
ble only for very low helium charge, much lower than derived 
from NLTE equations. We discussed the possibilities that may 
change this result, however, we were not able to find any con- 
sistent model that could lead to helium decoupling in the wind 
under such conditions. 

The main result of our paper, that for realistic values of he- 
lium charge helium will not decouple from the wind flow for 
parameters suitable for a star like cr Ori E, is in contradiction 
to the result of Hunger & Groote ( 1999 1 who obtained the op- 
posite behaviour. The reason for this is that helium is strongly 
coupled to hydrogen and this coupling cannot be neglected in 
model calculations. 

We have shown that for lower metallicities the wind mass- 
loss rate is significantly lower. This is caused by a signifi- 
cant decrease of the number of optically thick lines. For such 
low mass-loss rates the decoupling of wind components and 
a subsequent fall back of the passive component may occur. 
However, helium and hydrogen remain coupled in this case, 
consequently this does not explain the chemical peculiarity. 

We studied the multicomponent flow analytically and de- 
rived a simple condition for helium decoupling in the stellar 
wind. This analysis supports our basic result that helium decou- 
pling is possible only for a very low helium charge (i.e. when 
helium is basically neutral). 

We have also studied the possibility of decoupling of neu- 
tral helium. This could be an interesting way of helium de- 
coupling (at least partially), since the cross-section of colli- 
sions of protons with neutral helium is much lower than the ef- 
fective cross-section for collisions between charged particles. 
However, the neutral helium atoms do not remain neutral for 
sufficiently long time since the characteristic helium ionization 
time is much shorter than the characteristic time of change in 
flow properties Tfl OW - 
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As observations seem to favour the hypothesis of helium 
decoupling (Groote & Hunger 1997 Hunger & Groote 1999 1 
our result leaves us in the very unpleasant situation that we do 
not have any self-consistent explanation for the chemical pe- 
culiarity of He-strong stars. While the scenario proposed by 
Hunger & Groote ( 1999 1 requires artificial lowering of the he- 
lium charge, models for helium diffusion moderated by the stel- 
lar wind (Vauclair 1975 Michaud et al. 1987 1 require artificial 
lowering of the wind mass-loss rate. The only remaining influ- 
ence on the helium decoupling, which is very difficult to cal- 
culate is that of a strong magnetic field. Nevertheless, it seems 
clear that decoupling of helium in normal (non-magnetic) stars 
with comparable stellar parameters can be excluded as a mech- 
anism of creating the helium overabundance in the stellar atmo- 
sphere. Moreover, we were not able to find any self-consistent 
mechanism that would enable the decoupling in the case of 
magnetic stars. 

Appendix A: Approximate description of H-He-i 
flow 

In order to obtain an approximate description of a well-coupled 
multicomponent flow of hydrogen, helium and heavier ions 1 
(neglecting electrons) we can write the momentum equations 
of hydrogen (H) and helium (a) as (Burgers 1969 neglecting 
gas-pressure term and polarization electric field) 



V'h- 



dVH 

dr 



dr 



Kua „ , , Xu a Km XiH 

Gch(XHa)-. r + GchixiH): 

PH l*Horl Ph Rh| 

(A. la) 

K-Ha ^ , x X\i a K la X\ a 
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We assume that the flow is well-coupled, vh^ ~ v »77 = y 37 
and we use only the linear term of the Taylor expansion of the 
Chandrasekhar function Gch(x) 
Eqs. iA.l\ as 



^4. Then we can rewrite 
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(A.2b) 



where we used x a b - Av^/a^. Here k a b = 2K a b/ (3 -Jn\ and 
the velocity difference Av a b = V a - Vb for all a, b, which stand 
for H, a, and i. Subtracting Eqs. (IA.2b> and (IA.2a> we obtain 
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Due to the very low number density of heavier ions n u we can 
neglect the second right hand side term and derive (assuming 
that the temperatures of all wind components are equal; note 

1 Note that in this Appendix all heavier ions are described as one 
component. 



that the similar term on the left hand side can not be neglected 
because all terms on the left hand side are of the same order) 
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The sign of Av'H ff depends on the sign of the first bracket in 
Eq. dA~4l For 



i - mz j > 0, 



(A.5) 

L a lJ -ia<-Y^ 

i.e. for a low helium charge, the hydrogen velocity is higher 
than the helium velocity. Because n\ <sc n a we derive Av^ » 
Av'H ff for a well coupled flow and realistic values of the other 
parameters. Helium decoupling may occur in the case Avn a » 
AviH. Apparently, this may occur for very low helium charge 
z a , as was demonstrated in the paper. 

Finally, let us compare the frictional force acting on helium 
due to heavier ions and hydrogen in the case of a well-coupled 
flow. The frictional acceleration after Eq. (|3J, assuming > 
v H > V a and Avj H » Avhq., is 

1 1 

— K m Gca(Xi a ) H Kft a Gch(x\\a) ~ 

Pa Pa 



Av iH /, Av Ha ^ 



Pa a\a 



ha an a Avm , 



(A.6) 



Inserting Eq. iA.4\ for the fraction of velocity differences and 
performing just an order of magnitude estimate we may assume 
a xa ~ ana, Zi ~ z a and an a w am and conclude that both terms 
in 

kua a ia Av Ha 



1 + 



(A.7) 



ha au a Av'i H 

are of the same order. This means that the contribution of heav- 
ier ions and hydrogen to the helium acceleration are of the same 
order in the case of a well-coupled flow (e.g. for O supergiants). 
However, in the case when heavier ions are not able to effec- 
tively accelerate helium the friction with hydrogen dominates. 
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